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The superconducting (SC) state of iron-based compounds in both tetragonal and orthorhom- 
bic phases is studied on the basis of an effective Hamiltonian composed of the kinetic energy 
including the five Fe 3d-orbitals, the orthorhombic crystalline electric field (CEF) energy, and 
the two-orbital KugeF-Khomskii-type superexchange interaction. Our basic assumption is that 
the antiferromagnetic (AF) state in the parent compounds can be described by the d xz and 
d yz orbitals, and that the electrons in these orbitals have relatively strong electron correlation 
in the vicinity of the AF state. In order to study the physical origin of the structure-sensitive 
SC transition temperature, the effect of orthorhombic distortion is taken into account as the 
energy-splitting, A ort ho. , between the d xz and d yz orbitals. We find that the eigenvalue of the 
linearized gap equation decreases accompanied with the reduction of the partial density of 
states for the d xz and d yz orbitals as A ort ho. increases, and that the dominant pairing symme- 
try is an unconventional fully gapped s+_-wave pairing. We also find large anisotropy of the SC 
gap function in the orthorhombic phase. We propose that the CEF energy plays an important 
role in controlling T c and the SC gap function, and that orbital-selective superconductivity is 
a key feature in iron-based superconductors, which causes the structure-sensitive T c . 

KEYWORDS: multiorbital system, structure-sensitive superconducting transition temperature, iron-based 
superconductor, orbital-selective superconductivity, orthorhombic distortion 



1. Introduction 

Since the discovery of iron-based superconductors, 1 -* 
many works have been carried out on the physical prop- 
erties of these new superconductors both theoretically 
and experimentally 2 ) Except for in recent studies, 3 ' 4 ) 
superconductivity has been observed in the tetragonal 
phase, which is the crystalline structure above the struc- 
tural transition temperature T3. 2 ) The crystalline struc- 
ture becomes orthorhombic below Ts. Extensive research 
on these iron-based compounds has clarified several gen- 
eral features: a structure-sensitive superconducting (SC) 
transition temperature T c , 5 ' 6 ) the existence of discon- 
nected Fermi surfaces (FSs), 7 ~ n ) a small difference be- 
tween Tg and the Neel temperature Tn, 12 ~ 14 ) a small 
magnetic moment per Fe atom, 12,15 ' 16 ) a stripe-type 
antiferromagnetic (AF) spin structure, 12, 15 ' 17 ' 18 ) spin- 
singlet pairing, 18 ~ 22 ) various power-law dependences on 
the temperature of the spin-lattice relaxation rate T-f 
below T c , 18 - 29 ) the absence of the Hebel-Slichter (HS) 
peak just below T c , 18 ~ 29 ) and the robustness of T c against 
nonmagnetic impurities. 30 ~ 32 ) Note that the resistivity 
above Tn is metallic; 2 ) thus, it is appropriate to regard 
the magnetically ordered state a spin-density-wave state. 
For simplicity, we call this magnetically ordered state an 
AF state in the following. The phase diagrams of 1111, 
122, 11, and 111 systems are similar. The only difference 
is the following: in the 1111 system, the SC phase and 
AF phase are separated, 1 ) while the two phases adjoin 
each other in the other systems. 33 ~ 35 ) 

In early theoretical studies, 36,37 ) it was proposed that 



the symmetry of the SC gap function is an unconven- 
tional s-| -wave. Several experimental results have sup- 
ported this proposition, for example, the existence of a 
resonance peak in neutron scattering, 38 ) the absence of 
the HS peak of Tf 1 just below T c , 18 ~ 29 ) its revival in 
an overdoped compound, 39 ) and quasi-particle interfer- 
ence patterns. 40 ) On the other hand, experiments on the 
nonmagnetic impurity effect 30-32 ) indicate that the iron- 
based superconductors are robust against nonmagnetic 
impurities. This is inconsistent with the above theoret- 
ical predictions, since some theoretical studies on the 
impurity effect based on a standard T-matrix approx- 
imation using realistic parameters 41 ) have shown that 

T c with the s_| -wave symmetry decreases more rapidly 

than the experimental results. 30 ~ 32 ) Although it has been 
proposed 42,43 ) that the orbital fluctuation plays an im- 
portant role in explaining this discrepancy, the pairing 
mechanism of iron-based superconductors is still contro- 
versial. 

We consider that the key to clarifying the pairing 
mechanism is to understand the role of lattice distor- 
tion in determining the electronic states since the un- 
usual properties will be structure-sensitive. In this sense, 
the d xz and d yz orbitals will be more important than the 
other Fe 3d orbitals since they are closely connected to 
the orthorhombic distortion. The results of several re- 
cent experiments for the underdoped 122 system have 
suggested the importance of the orbital degree of free- 
dom and lattice distortion: anisotropic in-plane resistiv- 
ity above Ts, 44 ) anisotropic band dispersion above Ts, 45 ) 
the large incoherent electronic Raman spectrum for the 
d xz and d vz orbitals, 46 ) the smaller mobility of holes ob- 
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served in the Hall resistivity measurement using a two- 
carrier model, 47 - 1 the orbital-dependent modification of 
the electronic states across Tn observed in the angle- 
resolved photoemission spectroscopy (ARPES) measure- 
ment, 48 ) the marked softening of the elastic constants 
at Ts observed in ultrasound measurements, 49 ' 5 °) the 
good correlation between the orthorhombic distortion 
and the AF order observed in neutron and X-ray mea- 
surements. 34, 51 ) It is probable that the AF state in the 
parent compounds can be described by the d xz and d yz 
orbitals. 52 " 54 ) In both parent and underdoped 122 and 
11 systems, a large enhancement of the effective mass 
has been observed in several experiments. 8 " 11, 55 ~ 60 ) 

From the above experimental results, we expect that 
the electron correlation in the d xz and d yz orbitals will 
play an important role in both parent and underdoped 
compounds. Although there have been several theoretical 
proposals 52, 61 ~ 63 ) that a coupled spin and orbital order 
occurs in the AF state, the effect of this coupling on the 
SC state has not been clarified yet. 

In this paper, we focus on the SC state in both tetrag- 
onal and orthorhombic phases in order to investigate the 
effect of orthorhombic distortion on T c and the role of 
the spin-orbital coupling in determining the SC state. 
We assume that relatively strong electron correlation ex- 
ists for the d xz and d yz orbitals, and that the super- 
conductivity is induced from the Kugel'-Khomskii (KK)- 
type superexchange interaction 64 ) for the d xz and d yz 
orbitals. We call this orbital-selective superconductivity 
in the sense that only the d xz and d yz orbitals induce 
superconductivity. For this purpose, we introduce an ef- 
fective model composed of the kinetic energy including 
the five Fe 3(i-orbitals, the orthorhombic crystalline elec- 
tric field (CEF) energy, and the KK-type superexchange 
interaction. There have been some previous studies that 
discuss the SC state induced from the superexchange in- 
teraction. 65 ) In the present paper, a more realistic model 
Hamiltonian is considered. We use the two-dimensional 
(2-D) tight-binding model of the five Fe 3d orbitals down- 
folded from the local density approximation (LDA) cal- 
culation. 37 ) In order to study the effect of orthorhombic 
distortion, we introduce the splitting of energy levels, 
A or tho., between the d xz and d yz orbitals. Furthermore, 
we will include the effect of Coulomb interaction by a 
band renormalization in which the total bandwidth of 
the five Fe 3d-orbitals and the hybridizations concerning 
the d xz and d yz orbitals is reduced. This procedure is 
based on the same principle as the Gutzwiller approxima- 
tion. 66 ' From the experimental result 7 " 11,55 " 60,67 ) that 
the effective mass is at least twice the band mass, we use 
a renormalized band whose width is nearly half that of 
the original one. 

We study the behavior of the eigenvalue, A e , of the 
linearized gap equation and the pairing symmetry for 
< Aortho. < 0.1 eV. We find that A e decreases as 
A ort ho. increases, and that the SC state in the tetrag- 
onal phase gives the maximum T c . The former behavior 
is accompanied by the reduction of the partial density of 
states (pDOS) for the d xz and d yz orbitals. We also find 
that the dominant pairing symmetry is a fully gapped 



-wave pairing both in tetragonal and orthorhombic 

phases, and the second most dominant pairing symme- 
try is a d xy -wave pairing whose A c rapidly decreases as 
A or tho. increases. The finite hybridizations between the 
strongly correlated orbitals and the weakly correlated or- 
bitals cause the change in the pDOS for the d xz and d yz 
orbitals near the Fermi level. This is one of the physical 
origins of the sensitivity of T c . Large anisotropy of the 
SC gap function is found in the orthorhombic phase. 

The outline of this paper is as follows: In §2, we explain 
the form of the kinetic energy and the method used to 
include the effect of orthorhombic distortion, and we de- 
rive the KK-type superexchange interaction for the d xz 
and d yz orbitals. In order to discuss the SC state for the 
spin-singlet pairing, we use the mean field approximation 
(MFA) in our effective model and derive the linearized 
gap equation. §3 is devoted to showing the results of the 
eigenvalue and the pairing symmetry for various values 
of A ort ho. • In §4, we address the physical meaning of the 
obtained results, compare them with other previous the- 
ories, and discuss their correspondence with previous ex- 
perimental results. The paper concludes with a summary 
of our results in §5. 

2. Formalism 

In this section, we explain the methods used to take 
account of the effect of orthorhombic distortion and to 
construct an effective Hamiltonian, -ff ff, for the discus- 
sion of iron-based superconductors in both tetragonal 
and orthorhombic phases, assuming that the d xz and d yz 
orbitals have relatively strong electron correlation. The 
effective Hamiltonian is 

Hefi = H o + -f/ortho. + #int , (1) 

where Hq is the kinetic energy of the tetragonal phase 
modified by band renormalization, the detail of which is 
described in §2.1, H OItno , is the orthorhombic CEF en- 
ergy, and Hi n t is the effective interaction. In the follow- 
ing, the coordinates x and y are chosen in the directions 
of the nearest-neighbor (n.n.) Fe atoms in a unit cell. 
Note that, in the orthorhombic phase, the a-axis corre- 
sponds to the ^-direction. For convenience, the five Fe 
3<i-orbitals, d xz , d yz , d xy , d x 2_ y 2, and d z 2, are labeled 
1, 2, 3, 4, and 5, respectively. We define the band fill- 
ing n e as the number of electrons per site. In iron-based 
compounds, n c is 6 + x, where x represents the doping 
level. 

2.1 Kinetic energy and orthorhombic CEF energy 

In order to describe the electronic states of iron-based 
superconductors, we use the 2-D five-orbital model down- 
folded from the LDA calculation, 37 ) 

h«= E E E(4 j 4 CT ^ + h.c.) 

i.j a, 6=1 a 

5 

+ E E E(^a - H)n ia a , 
i a—1 a 

5 

= E E Y.tab{k)cl aa c kba , (2) 

k a, 6=1 a 
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(a) Tetragonal phase 



(b) Orthorhombic phase 
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Fig. 1. Schematic energy levels for the five Fe 3d orbitals in the 
(a) tetragonal and (b) orthorhombic phases. 



where c\ acr (c iaa ) is the creation (annihilation) operator 
that creates (annihilates) an electron in orbital a with 
spin a at site i, h.c. means the Hermitian conjugate, 
niaa = c\ aa c iaa , c kaa is the Fourier component of c iaa , 
and t* j , E a , and e a b(k) denote the values of in-plane hop- 
ping integrals, on-site energies, and energy dispersions 
measured from the chemical potential fi, respectively. 

The effect of electron correlation is taken into account 
by modifying H to 



H o = E E H^ab{k)cl aa c kbo 

k a, 6—1 u 



(3) 



where e a b(k) is the renormalized energy dispersion. 
e ab (fc) is constructed as follows. All t*£ in eq. (2) are re- 
duced by a factor of 0.6, assuming band renormalization 
due to the electron correlation effect. In addition to this, 
t\'£ and tgo ( a = 1 — 5) are further reduced by a factor 
of 0.75. We use this modification assuming a relatively 
strong electron correlation in the d x 
Finally, the on-site energies, E a (a =1 
as 



. , z and d yz orbitals. 
5), are adjusted 



E a 



E a + 0.0045 (eV) for a = 1, 2, 3 
E a - 0.03 (eV) for a = 4, 5 



(4) 



in order to make the total occupation number of the d xz 
and d yz orbitals, {n c ) xz+ y Z , close to 3, since we assume 
a KK-type interaction as discussed below. Note that the 
center of the energy levels remains zero under the trans- 
formation in cq. (4). 

In order to include the effect of orthorhombic distor- 
tion, we use the orthorhombic CEF energy for the d xz 
and dyz orbitals (i.e., a = 1 and 2), 



^ortho. = 2-, E o \ n k2o 



n k ic 



(5) 



where A ort ho. is the splitting between the energy levels 
of the d xz and d yz orbitals. It has been experimentally 
shown that the a-axis is longer than the &-axis in the 
orthorhombic phase. 2 ) Since As ions have negative charge 
and the d xz orbital extends toward As ions, we expect 
that the energy level of the d xz orbital is lower than that 
of the dyz, orbital, i.e., 

E X ^E X - E 2 ^E 2 + (6) 

Note that this expectation is consistent with the result 
of an ARPES measurement. 45 ' Schematic energy levels 
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Fig. 2. Schematic diagrams of possible intermediate states on a 
single site and the corresponding energies neglecting the effect of 
Aortho. - The upper and lower horizontal lines represent the d yz 
and d xz orbitals, respectively. Up (down) arrows correspond to 
spin-up (spin-down) electrons. 



for the five Fe 3d orbitals are shown in Fig. 1. 

In this work, we use the value of A ort ho. as a parameter 
ranging from to 0.1 eV. According to the ARPES mea- 
surement, 45 ' A ort ho. is approximately equal to 0.06 eV at 
T/Xn ~ 0.58 for BaFe 2 As 2 . Moreover, the LDA calcula- 
tion for the parent compounds 68 -* indicates that A ortrio . is 
on the order of 0.1 eV. In order to maintain the topology 
of the FS, we choose A ortri0 . so as to satisfy the inequali- 
ties A ortho . < 2(E 1 - E 5 ) and A ortho . < 2(E 3 - E 2 ). Pre- 
cisely speaking, the orthorhombic distortion also leads 
to anisotropies in the hopping integrals. However, we ne- 
glect this effect in the following calculations, the validity 
of which will be discussed in §4. 

Ho + Aortho. is diagonalized as 



H + #ortho. = E E E l a{k)c\ aa C ko 
k a a 



(7) 



where a is the band index and e a (k) is the renormalized 
energy dispersion for band a. Here, c kaa is related to 
Cfcacr by the unitary transformation 



aa c kc 



(8) 



2.2 KK-type superexchange interaction 

We derive the KK-type superexchange interaction be- 
tween the d xz and d yz orbitals in order to discuss the 
electronic states of iron-based superconductors. We as- 
sume that relatively strong electron correlation exists in 
the d xz and d yz orbitals since the AF state is realized 
mainly in these orbitals. 52 ~ 54 ) 

First, we assume standard on-site interactions for the 
d xz and d yz orbitals, i.e., the intraorbital Coulomb inter- 
action U, the intcrorbital Coulomb interaction U', and 
the Hund's rule coupling J H . (The pair-hopping term 
does not directly affect the electronic state for the case 
with (n e ) xz + yz = 1 or 3 in the strong coupling limit.) In 
this work, U, U' > Jh is assumed. 

Assuming that (n c ) xz+yz is close to 3 or 1 and us- 
ing the second-order perturbation theory from the strong 
coupling limit, we obtain 



2 

Hint = E E E 

i,u j— n.n. a, a', 6—1 



,t J 



C ia'a C jbs\^J 



(CjbaCiaa C-jboC-iaa) 



I v-> tabtb'a' f f , _ n 

' 2-i 2-i 2-i nrr, C ia'cr C jb>g-^ C jba c iaS CjbgC-iacr) 

i,cr j=n.n. a,a',b,b' = l ^ u 

. v v v ( -ni+fr+fr' tabtb ' a ' J J_ 

i,<7i=n.n.o,o',6,6' = l 2 \ U ~ J K) 3 

( c jba C iaa + c j oS Ci acr ) 
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+ E E E (-i) 1+b+b ' 

i,a- j—n.n. a,a',b,b f — 1 

+ (n.n.n. terms), 



tabtb'a' ft- 

T~ C ia'a C jb< cr C jba c iaa 



u 



(9) 



where a and cr represent orbital and spin degrees of free- 
dom, respectively, j =n.n. denotes the summation over 
the n.n. sites of i, a represents the other orbital from 
a, and a — —a. Possible intermediate states and corre- 
sponding energies are shown in Fig. 2. In eq. (9), we have 
also included the superexchange interactions originating 
from the next-nearest-neighbor (n.n.n.) hoppings, since 
they are not negligible in iron-based compounds. Actu- 
ally, the magnitudes of the hopping integrals are strongly 
affected by hybridizations between the 2p orbitals of As 
(or P, Te, or Se) and the 3d orbitals of Fe. 

When (n c ) xz+yz is away from 3 or 1, many other su- 
perexchange interaction terms appear. In the doped case, 
these terms may affect the properties of the SC state. 
However, we neglect these effects as far as (n c ) xz+yz is 
close to 3 or 1. Furthermore, we have neglected the effect 
of A ort ho. in the denominators in eq. (9) since we assume 
A or tho. << U, U', (V — Jh), with the latter being on the 
order of 1 eV. Note that we do not use the pseudospin 
and spin operators in the KK-type interactions since we 
intend to study the SC states in the doped cases. 

We consider that the KK-type superexchange inter- 
action in eq. (9) can describe the stripe-type AF state 
when there is no doping or when (n c ) xz+yz — 3. When 
an orbital-ferromagnetic phase transition occurs, our ef- 
fective interaction becomes similar to the anisotropic 
Hcisenberg spin Hamiltonian with n.n. and n.n.n. inter- 
actions. 69 ' 70 - ) In this case, the stripe-type AF state can be 
described naturally. Furthermore, the tight-binding cal- 
culation for the 1111, 122, and 111 systems shows that 
the magnetic interactions are short-range. 71 ' This is con- 
sistent with our KK-type interaction. 

2.3 Superconductivity 

As discussed in §1, we assume that the KK-type su- 
perexchange interaction for the d xz and d yz orbitals in 
eq. (9) induces superconductivity, and thus apply the 
MEA to these interaction terms. Since the terms with a 
denominator of (V — Jh) in eq. (9) give no contribution 
to the spin-singlet pairing, we consider only the terms 
with U and U' in the denominator. Note that there are 
hybridizations between the five Fe 3<i-orbitals, while the 
superexchange interactions are only between the d xz and 
d yz orbitals. 

After a straightforward calculation, the effective 
Hamiltonian within the MEA is obtained as 

^cT A = E E £Mfc)cjL<*fc, 

k a, 6—1 o 
. ^ortho. / \ 

+ E E — o — ( nfc2 <T -n k ia) 

k cr * 

+ E E (Aat(fc) C fcat C -febi + h.C.). (10) 

k a,b=l 

Here, A a b(k) is the SC gap function for the d xz and d yz 



orbitals defined as 
1 



A ab (k) 



N 



E E v 

k' a'.b'- 



(ii) 



with N being the number of lattice sites of Fe atoms. 
For simplicity of the numerical calculation, we introduce 
the cutoff energy e c , i.e., the fc-sum is restricted. Here, 
V aa ',bb'{k — fe') is the coefficient of the two-body interac- 
tion between (fe a t, —kb\) electrons and (fe' a' t, —k' 
b' I) electrons, which is obtained from eq. (9) by Fourier 
transformation. Explicitly, V a ' a " ,vb" (fe - k') is written as 

W, y&»(fc - fe') = Ja'a",Vb»f(k - fe'), (12) 

where f(K = k — fe') and J a 'a",b'b" are defined as 

^cos^) for (i-j) = (1,0) and (-1,0) 
f(K) = I 2cos(K y ) for (i-j) = (0,1) and (0,-1) , 
4cos(K x ) cos (K y ) for (i — j) = n.n.n. 

(13) 

and 



Ja'a" ,6' i 



At 2 

(*ii+*2 2 ) ~ 4t 2 



U' 



zt 12 

U' 

+ 



2t 2 



12 

u 



12 



= < 



2t U 't 

^ I ll I 22 

w 

2t aa tl2 ^12(^11+^22) c I, 

+ w forfl 

2t aa t\i t\i(t\\ + £22) , 1 

+ W fOTfl 

2t aa t\i £12(^11 +£22) , 1 
+ W f ° ra 

2t aa tl2 , ^12(^11+^22) r 1 

H tt. for a 



U 



U' 



for a' = b' = a" = b" = a 

for a' = a" = a,b' = b" = a 

for a' = b' = a, a" = b" = a 

for a' = b" = a,b' = a" = a 

b" = b' = a, a' = a 

a" = b' = a, b" = a 

a" = b" = a,b' =a 

b" = b' = a, a" = a 
(14) 



respectively. Note that the hopping integrals in, i22, and 
ti2{= ^21) depend on the direction of (i — j). The form 
of J a 'a",b'b" indicates that the number of SC attractive 
channels for electrons is enhanced by the coupling be- 
tween spin and orbital degrees of freedom. 

At T — T c , the SC gap function satisfies the following 
self-consistent equation: 



A a6 (fc) 



1 

N 



E 



2 

E 



k' o',6',o",6" = l a' 

eV(fe') , tanh !^) ); 



E v aa ,, bb ,(k-k') 



&a»b»(k') 



(tanh 2T C ' — 2T C >2{e a ,(k') + i,A-k>)) 

a'a"(U-k')b'P'(U}_ k ,)p>b"- (15) 

Note that the sum of the orbital indices (a', b' , a", b") is 
restricted to the d xz and d yz orbitals, while there is no 
restriction for the sum of the band indices (a',/3'). The 
effect of the other three Fe 3<i-orbitals (i.e., d xy , d x 2_ y 2, 
and d z 2) is included in the energy dispersions and uni- 
tary matrices, since the d xz and d yz orbitals have finite 
hybridizations with the other three Fe 3d-orbitals. T c can 
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be estimated by the following linearized gap equation: 

E v aa , M ,(k-k') 

ot>,P< 



A c A ab (fc) = i-E E 

iv fc' a',6',a",fc" = l a'jjS' 



( tanh ■ 



e /3' 



-fc')) 

(16) 



2T ' 2T '2(60,/ (fc') 

(DiiOa'«'(^)a'a»(^-fcO^'(^ifc/W. 
where T c corresponds to the temperature at which the 
eigenvalue A becomes unity. 

We remark on the symmetry of the pair amplitude in 
terms of the orbital degree of freedom. The pair ampli- 
tude satisfies 



(c-kaic kbt ) = ±(c_ fe6J .Cfe ot ) for a ^ b , 



(17) 



where the upper and lower signs correspond to orbital- 
ferromagnetic and orbital- AF pairings, respectively. Ac- 
cording to the form of the interaction given by eq. (14), 
the attractive interaction is dominant in the case of 
orbital-ferromagnetic pairing. In the following calcula- 
tions, we thus consider only the orbital-ferromagnetic 
spin-singlet pairing. 

3. Results 

In this section, we show the results of numerical cal- 
culations using the MFA. We take 254x254 meshes in 
the unfolded Brillouin zone (BZ) and set n c =6.24, 
e c =0.01 eV, (U,U',J U ) = (1,0.8,0.1) (eV), and the 
temperature T =0.003 eV. For convenience, the SC gap 
function is always normalized. According to the 2-D 
tight-binding calculation in the tetragonal phase, 37 ) hop- 
ping integrals for the d xz and d yz orbitals are given by 
(*?i n -)||x = m?-)\\ v = 0.08 eV, (V£t) h = (t^-)\\ x = 



0.34 eV, i? 2 n n - 
and i? 2 n n - = 



O.OOeV, t 



li 



j-ii.n.n. 

— ^22 



= -0.24 eV, 



0.09 eV. (e.g., (1%t)\\x denotes (i*f ) for 
(i — j) = (±1,0), etc.) Hereafter, the unit of energy is 
taken as 1 eV. 

We remark on the validity of our choice of Jh and 
U. (V is given by the relation U' = U — 2Jh.) A recent 
study on the five-orbital Hubbard model within the MFA 
shows that the small magnetic moment in the parent 
compounds can be understood only in the case of small 
values of Jh/U. 72 ^ This is also supported by Gutzwiller 
approximation studies on the two- and three-orbital Hub- 
bard models. 54,73 ^ In accordance with these results, we 
use small values of Jn/U to discuss the SC state in the 
vicinity of the AF state. 

3.1 Renormalized band structure 

As described in §2.1, a renormalized band structure is 
constructed in order to take into account the relatively 
strong electron correlation for the d xz and d yz orbitals. 
By this procedure, (n c ) xz+yz becomes larger than that 
in the original band structure and is equal to 2.62, which 
is close to 3 (Table I). 

The obtained band structure, the corresponding FSs, 
and the density of states (DOS) in the tetragonal phase 
are shown in Figs. 3-5, respectively. For comparison, the 
original band structure, the corresponding FSs, and the 
DOS are also depicted in the figures. 



Table I. Occupation numbers of electrons of the five Fe 3d or- 
bitals in the original and renormalized band structures in the 
tetragonal phase. 



Orbital 


Original 


Renormalized 


dxz 


1.19 


1.31 




1.19 


1.31 


d X y 


1.14 


1.20 


dx 2 —y 2 


0.98 


0.98 


d z 2 


1.49 


1.46 



(a) 



(b) 




(0,0) 



Fig. 3. (Color online) (a) Original band structure and (b) renor- 
malized band structure in the tetragonal phase. The dashed lines 
represent the chemical potentials. 
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Fig. 4. FSs for the (a) original and (b) renormalized band struc- 
tures in the tetragonal phase. 
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Fig. 5. (Color online) DOSs for the (a) original and (b) renormal- 
ized band structures in the tetragonal phase. The dashed lines 
represent the chemical potentials. 



Similarly, we calculate the FSs and DOS in the or- 
thorhombic phase. As a typical case, the results for 
A or tho. = 0.08 arc shown in Fig. 6. The main effects of 
orthorhombic distortion are i) the disappearance of one 
of the FSs around the T-point [i.e., k = (0,0)], and ii) 
the reduction of the DOS in the vicinity of the Fermi 
level. 

3.2 Eigenvalue of the linearized gap equation 

We investigate the eigenvalue, A c , of the linearized gap 
equation, eq. (16), for < A ort ho. < 0.1. Figure 7 shows 
the obtained eigenvalue, A e , as a function of A ort ho .• We 
find that A e decreases as A ort ho. increases. This result 
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Fig. 6. (Color online) (a) FS and (b) DOS for A ortho . 
The dashed line represents the chemical potential. 
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Fig. 10. (Color online) Band structures around the T point for 
^ortho. = ( a ) 0.02, and (b) 0.04. Solid lines and dashed lines 
correspond to the chemical potentials and cutoff energies e c (= 
0.01), respectively. 
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Fig. 7. (Color online) Eigenvalue of the linearized gap equation, 
A c , for < A ortho . < 0.1. 



Fig. 11. (Color online) Eigenvalue of the linearized gap equation 
A e with s-wave and d-wave symmetries for < A ort ho. < 0.1. 
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Fig. 8. (Color online) (a) NOS and (b) pNOS for d xz and d yz 
orbitals in the vicinity of the Fermi level as a function of A ortrio . 
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Fig. 9. FSs for A ortho . = (a) 0.02, and (b) 0.04. The main changes 
in the FSs are shown by arrows. 



means that T c is highest in the tetragonal phase. In order 
to investigate the physical origin of this behavior of A e , 
we calculate the number of states (NOS) near the Fermi 
level, which is defined as the DOS integrated from — e c to 
£ c . Figure 8 shows the total NOS and the partial number 
of states (pNOS) for the d xz and d yz orbitals. Comparing 
Fig. 7 with Fig. 8, we consider that the physical origin of 
the behavior of A c is the variation of the pDOS for the 
d xz and d yz orbitals. This result indicates that energy- 
splitting due to the orthorhombic distortion plays an im- 
portant role in controlling T c . A similar conclusion was 



(a) A n 
jt 
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Fig. 12. (Color online) SC gap functions An and A22 for 
Aortho. = 0. The dashed lines represent the FSs for A ortrio = 0. 



also obtained theoretically for CeCoIns. 74 ) 

Figure 7 also shows that A e decreases with increas- 
ing A ort ho. for small values of A ort ho. , while it does not 
change greatly for A ortho . > 0.05. We speculate that this 
difference originates from the existence of a small hole 
pocket around the T-point. Figures 9 and 10 show the 
main changes in the FSs for A ort ho. = 0.02 and 0.04 and 
the corresponding band structures near the T-point, re- 
spectively. The arrows in Fig. 9 represent the changes 
in the FSs. One of the small hole pockets around the 
T-point seems to disappear near A ort ho. = 0.04. We con- 
sider that the small hole pocket enhances A . When it 
vanishes above A or tho. = 0.04, A e does not change fur- 
ther. A similar effect has been observed in a previous 
study for a system with multiple FSs. 75 ) 

3.3 SC gap function and pairing symmetry 

We now turn to the SC gap function and its pair- 
ing symmetry. Before we show the results of numerical 
calculations, we remark on the method used to deter- 
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Fig. 13. (Color online) SC gap functions An and A22 for 
Aortho. = 0.04. The dashed lines represent the FSs for A or tiKj. = 
0.04. 
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Fig. 15. (Color online) (a) Maximum values of orbital-diagonal 
SC gap functions for s+_-wave symmetry, max An and max 
A22, as a function of A ort ho. and (b) pNOSs for the d xz and d yz 
orbitals normalized by those in the tetragonal phase. 



mine the pairing symmetry for the multiorbital SC state 
within the linearized gap equation. 74 ) For the multior- 
bital SC state, there are two components of the SC gap 
function: orbital-diagonal and orbital-off-diagonal com- 
ponents. The pairing symmetry is determined from the 
orbital-diagonal SC gap function since it belongs to the 
same irreducible representation of the SC state. Here- 
after, we show only the orbital-diagonal SC gap func- 
tions. 

First, we investigate the most dominant pairing sym- 
metry. As shown in Fig. 11, an unconventional fully 

gapped S-| -wave pairing is the most stable for < 

A or tho. < 0.1. As typical forms of the SC gap functions, 
An and A22 in the tetragonal and orthorhombic phases 
are shown in Figs. 12 and 13, respectively. At A ortho = 0, 
the second most stable pairing symmetry is a d^y-wave 
pairing whose form is depicted in Fig. 14. Note that the 



values of A e for other d-wave pairings (e.g., the d x 2_ y 2- 
wave pairing) are smaller than that for the dry-wave pair- 
ing. 

Let us discuss the condition required to stabilize the 
d^-wave pairing. At A or tho. = 0, the FS around the X- 
or Y-point [i.e., (71-, 0) or (0, w)] does not have components 
of the d xz and d yz orbitals along the k x - or Ay-direction, 
respectively. 75 ^ Therefore, the energy cost is small even 
if a line node appears along the k x - and fc y -axes. This 
is why the d xy -wave pairing is stabilized. However, the 
rapid drop in A e for the d xy -wave pairing shown in Fig. 
11 may be due to a subtle balance of the FSs around the 
T-point. 

Finally, we investigate the anisotropy of the SC gap 
function in the orthorhombic phase. Figure 15(a) shows 
the maximum values of the SC gap functions, An and 
A22 , as a function of A ortri o. ■ When A or tho. is on the or- 
der of 10% of the renormalized bandwidth (A or tho. =0.1), 
the ratio of maxAn to maxA2 2 becomes 1:1000. Large 
anisotropy is thus expected for the SC state in the or- 
thorhombic phase. With increasing A or tho.i the ampli- 
tude of An decreases rapidly, while that of A 2 2 increases. 

[The SC gap functions for the -wave pairing become 

maximum at k = (0,0).] This behavior can be under- 
stood from the pNOS for the d xz and d yz orbitals shown 
in Fig. 15(b). We can see that A ort ho. markedly affects 
the pNOS for the d xz orbital, while it does not greatly 
affect that for the d yz orbital. This is because the energy 
level of the d xz orbital decreases to below the Fermi level. 
This variation of pNOSs leads to the decrease in An and 
increase of A 2 2 as A ortrio . increases. 

Since the experimentally observed SC gap function is 
the sum of the orbital-diagonal and orbital-off-diagonal 
components, the SC gap function becomes anisotropic in 
the case of different values of An and A22. Figure 15 
indicates that the SC gap function around the X-point 
becomes larger than that around the Y-point for the SC 
state in the orthorhombic phase, since A22 is dominant in 
the orthorhombic phase and has this anisotropy (Fig. 13). 
We thus expect a large anisotropy of the SC gap func- 
tion in the orthorhombic phase, which can be detected by 
ARPES measurement of the SC state in the orthorhom- 
bic phase. As we shall discuss in §4, this anisotropy is 
also expected in a SC state coexisting with a nematic 
state. 

4. Discussion 

First, we discuss the validity of neglecting the 
anisotropy of hopping integrals due to the orthorhom- 
bic distortion. According to a recent tight-binding cal- 
culation for the parent compound of LaFeAsO, 68 ) there 
is little significant difference in the hopping integrals 
for the d xz and d yz orbitals. Since we focus only on 
the region where the orthorhombic distortion is small 
(A or tho. < 0.1), the effect of the anisotropy on the hop- 
ping integrals is negligible. In other words, the effect of 
orthorhombic distortion on the energy-splitting will be 
more important than that on the hopping integrals in 
the presence of a small orthorhombic distortion. It is rea- 
sonable to consider only a small orthorhombic distortion 
since we consider the superconductivity in the doped case 
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where the effect of orthorhombic distortion is suppressed. 

We address the role of the spin-orbital coupling in de- 
termining the SC state. As shown in §2, the degeneracy 
of orbitals gives additional degrees of freedom which en- 
hance the number of attractive channels for Cooper pair- 
ing. On the other hand, the orbital degree of freedom 
also has a depairing effect on superconductivity. In the 
present MFA, we have not clarified the depairing effect. 

In iron-based compounds, As ions are located asym- 
metrically with respect to the Fe plane. 2 - 1 Therefore, 
there are many hopping integrals between Fe 3d orbitals, 
which are not allowed without As ions. This increases the 
number of interactions in eq. (9), which eventually in- 
creases the number of SC attractive channels. Note that 
these terms also enhance the magnetic instability; thus, 
they do not always favor the SC state. 

Next, we compare our theory with other previous the- 
oretical studies. As shown in §1, there are two differ- 
ent candidates for the pairing symmetry of the iron- 
based superconductors, i.e., s H -pairing 36 ' 37, 76 ~ 81 ' and 

s ++ -pairing. 42,43 ' Our theory suggests the former pairing 
symmetry. Although the mechanism is not clear yet, we 
consider that the antiferromagnetism is closely related 
to the emergence of the superconductivity Note that 
the quantum critical point of the antiferromagnetism has 
been observed for optimally doped BaFe2(Asi_ x P x ) 2 by 
a nuclear magnetic resonance (NMR) measurement. 82 ' 
This result indicates that the antiferromagnetism plays 
an important role in the emergence of superconductivity 
for iron-based compounds. 

We remark on the main difference between our the- 
ory and spin-fluctuation-mediated pairing, whether the 
glue of the superconductivity exists or not. In this pa- 
per, the SC attractive force originates directly from 
the KK-type superexchange interactions, while in spin- 
fluctuation theory, it originates from the fluctuation us- 
ing a wide range of frequency space. Our obtained results 
have only n.n. and n.n.n. SC gap functions. This is qual- 
itatively consistent with the SC gap functions obtained 
in the spin-fluctuation-mediated pairing, since the lat- 
ter can be reproduced by the n.n. and n.n.n. parameters 
in real space. 83 ' The SC gap function observed in sev- 
eral ARPES measurements can also be fitted with the 
n.n. and n.n.n. parameters in real space. 10,84 ' Therefore, 
the behaviors of A c and the anisotropy between An and 
A22 as a function of A ort ho. will be qualitatively general 
in iron-based superconductors. 

Next, we discuss the experimental correspondence for 
our assumption that the attractive interaction appears 
only among the special orbitals such as d xz and d yz . 
This assumption is based on the idea that the AF state 
in parent compounds can be described by two local- 
ized orbitals. 52 ' Then, electrons in these orbitals in- 
duce superconductivity due to chemical doping and/or 
physical doping. We consider that this orbital-selective 
AF state does not contradict the experimental results 
for the metallic behavior, since the other orbitals have 
itinerant character. (A similar situation is expected in 
Ca 2 -a;Sr 2 ;Ru04. 85 ') As shown in §1, there are several 
experimental results that support our assumption: the 
large incoherence of the electronic Raman spectrum for 



these orbitals, 46 ' and the smaller mobility of holes ob- 
served in Hall resistivity measurement using a two-carrier 
model. 47 ' Note that the hole mobility will mainly origi- 
nate from the d xz and d yz orbitals, since the hole pock- 
ets of the FS usually consists of these orbitals. More- 
over, theoretical calculations for FeSe using the dynami- 
cal mean field theory 86 ' 87 ' show that the renormalization 
factors for the d xz and d yz orbitals are small (= 0.28) 
compared with those of d z 2 and d x 2_ y 2. (Note, however, 
that the renormalization factor for the d xy orbital is also 
small.) This is consistent with our assumption. 

In Bai_ a; K :I ;Fe2As2 and BaFe2(Asi_ 2; P ; j;)2, however, 
the ARPES measurement 88 ' indicates that the d z i or- 
bital contributes to the FS, and the SC gap of the d z 2 
orbital is as large as that of the other orbitals. In this 
case, it is necessary to include the d z 2 orbital, although 
the d z 2 orbital can be a passive orbital for superconduc- 
tivity. 

We next address the possibility of superconductivity in 
the orthorhombic phase. Recently, the emergence of su- 
perconductivity in the orthorhombic phase with nearly 
100% volume fraction has been reported for single crys- 
tals of Bao.6Ko.4Fe2As2, 3 ' although there have been some 
reports for the same compound that claim the inho- 
mogeneous coexistence of SC and AF states. 26,89 ' Ac- 
cording to rcf. 3, T c increases as the orthorhombicity 
5 = (a — b)/(a + b) decreases. This is consistent with 
our results. As discussed in §3.2, the decrease in T c due 
to the small orthorhombic distortion can be understood 
as the reduction of the pDOS for the d xz and d yz or- 
bitals. Furthermore, in FeSe .g4 under high pressure, the 
microscopic coexistence of SC and AF states in the re- 
gion of higher pressure was observed by the muon-spin 
rotation measurement. 4 ' The SC state may also exist in 
the orthorhombic phase of Ba(Fei_ a; Co a ;)2As2. 33 ' 90 ' 

On the other hand, in SrFe 2 As 2 under high pressure, 
the NMR spectrum for the SC-dominant phase is almost 
equal to that for the tetragonal paramagnetic state. 18 ' 
This indicates that superconductivity is not realized in 
the orthorhombic phase. Although more detailed analysis 
in the presence of the AF order is needed, we consider 
that our effective model can qualitatively describe the 
properties of the SC state if superconductivity is realized 
in the orthorhombic phase. 

Finally, we remark on the nematic state observed 
in the 122 system, 44,45 ' since the SC state coexisting 
with the nematic state is nearly the same as that in 
the orthorhombic phase. The nematic state is charac- 
terized by the lowering of rotational symmetry without 
any changes in the crystal structure. There have been 
some proposals of the emergence of the nematic state 
in strongly correlated systems such as cuprate 91 ' and 
Sr 3 Ru 2 7 . 92 ' In the iron-based compounds, its emer- 
gence has been proposed on the basis of experimental 
results for Ba(Fei_ a; Co a ;)2As2 44 ' 45 ' and proposed theo- 
retically. 93 ~ 96 ' Although the effect of nematicity on the 
SC state has not been clarified, we propose that the 
anisotropy of the SC gap function, which is similar to 
that shown in Fig. 15(a), i.e., a difference between the 
SC gap functions around the X- and Y-points, can be 
observed if the SC state coexists with the nematic state. 
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5. Summary 

We investigated the eigenvalue, A , of the linearized 
gap equation and the pairing symmetry in a model for 
iron-based superconductors in both tetragonal and or- 
thorhombic phases on the basis of the MFA. Our effec- 
tive model consists of the five-orbital kinetic energy, the 
orthorhombic CEF energy, and the two-orbital KK-type 
superexchange interaction. 64 ) The effect of orthorhom- 
bic distortion on electronic states was taken into consid- 
eration as the energy-splitting of d xz and d yz orbitals, 
A or tho.- In order to take account of the relatively strong 
electron correlation, the procedure of band renormaliza- 
tion was used. We found a decrease in A c accompanied 
with a reduction of the pNOS of the d xz and d yz orbitals 
near the Fermi level as A ort ho. increases. This result is 
consistent with the results of a recent experiment. 3 ) The 

fully gapped s_| -wave symmetry pairing is predominant 

in both tetragonal and orthorhombic phases. A c for the 
d xy -wave pairing, which is the second most dominant 
symmetry, rapidly decreases as A or tho. increases. The 
highest T c for the fully gapped -wave pairing is ob- 
tained in the tetragonal phase. The SC gap function for 

the fully gapped -wave pairing becomes anisotropic 

in the orthorhombic phase due to the crystal symme- 
try. We found large anisotropy of the SC gap function in 
the orthorhombic phase, maxAn:maxA22 ^1:1000, al- 
though A ort ho. is on the order on 10% of the renormalized 
bandwidth. Anisotropy originating from the difference in 
the values of the dominant orbital-diagonal SC gap func- 
tion (i.e., A22) around the X- and Y-points can be de- 
tected by ARPES measurement if the superconductivity 
appears in the orthorhombic phase or in the coexisting 
phase with the nematic state. 

We propose that the CEF energy plays an important 
role in controlling T c and the SC gap function for mul- 
tiorbital superconductors, and consider that the orbital- 
selective superconductivity discussed in this paper will 
be universal for the structure-sensitive superconductiv- 
ity observed in multiorbital systems. 
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